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We study the weak link between current-carrying super- 
conductors, both conventional and d-wave. The state of the 
system is controlled by two parameters: the order parame- 
ter phase difference and the superfluid velocity Vs, which 
parameterizes the parallel to the boundary transport super- 
current which is injected externally. The low-temperature 
current-phase relations are derived. We consider two models 
of weak links: a constriction between two conventional su- 
perconductors and a plane boundary between two differently 
orientated d-wave superconductors. We show that for some 
relation between and Vs quasiparticles create the current 
along the boundary which flows in the direction opposite to 
the transport supercurrent. 



link. These quasiparticles create current in the direction 
opposite to the transport supercurrent in the banks. The 
aim of the present paper is to study this countercurrent 
theoretically. 

First we consider a weak link between two conventional 
superconductors as a constriction (see also [7]), which 
is the simplest situation to demonstrate how the coun- 
tercurrent appears and how the interface-induced quasi- 
particle states influence the current through the contact 
(Josephson current). We then consider a weak link be- 
tween two differently orientated rf-wave superconductors 
(see also [8]). The influence of the transport supercurrent 
on the state of such a system is interesting in the light of 
possible applications [9], [10]. 



I. INTRODUCTION 

It was demonstrated both experimentally [1]- [3] and 
theoretically [4] , [5] that in some high-Tc superconductors 
the paramagnetic Meissner effect takes place. Namely, 
at the boundary of a d-wave superconductor placed in 
external magnetic field the current flows in the direc- 
tion opposite to the diamagnetic Meissner supercurrent 
which screens the external magnetic field. This coun- 
tercurrent is carried by the surface-induced quasiparti- 
cle states. These nonthermal quasiparticles appear be- 
cause of the sign change of the order parameter along 
the reflected quasiparticle trajectory. Such a depairing 
mechanism is absent in the homogeneous situation. At 
T = in a homogeneous conventional superconductor the 
quasiparticles appear only when the Landau criterion is 
violated, at Vs > Aq/pf- (Here Vs is the superfluid ve- 
locity which parameterizes the current-carrying state, Aq 
stands for the bulk order parameter, and is the Fermi 
momentum.) 

In a weak link of two superconductors with order pa- 
rameter phase difference (j), the order parameter can be 
tuned by (j). In the constriction between two conven- 
tional superconductors the proximity order parameter at 
the contact is given by [6] A^ = Aq |cos((^/2)|. Suppose 
there is supercurrent tangential to the boundary between 
two clean conventional superconductors. Namely, let us 
consider a weak link between current-carrying supercon- 
ductors. For some values of phase difference (f) and super- 
fluid velocity when the condition VsPf > A^ is satis- 
fied, the quasiparticles appear in the vicinity of the weak 



II. THE MODEL. GENERAL REMARKS 

We consider a perfect contact between two clean singlet 
superconductors. The external order parameter phase 
difference (j) is assumed to drop at the contact plane at 
X ^ 0. The homogeneous supercurrent flows in the banks 
of the contact along the y-axis, parallel to the boundary. 
The sample is assumed to be smaller than the London 
penetration depth so that the externally injected trans- 
port supercurrent can indeed be treated as homogeneous 
far from the weak link. The size of the weak link is as- 
sumed to be smaller than the coherence length. Such 
a system can be quantitatively described by the Eilen- 
berger equation. Taking transport supercurrent into ac- 
count leads to the Doppler shift of the energy variable 
by p pVs ■ The standard procedure of matching the solu- 
tions of the bulk Eilenberger equations at the boundary 
gives the Matsubara Green's function G^{0) at the con- 
tact at X = [8]. The component Gl^{0) = g^{0) of 
Gi^{0) defines current density at the boundary: 



j(0) - AneNovpT (vlm.gc.(0))^ 
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where A^o is the density of states at the Fermi level, {■■■)^ 
denotes averaging over the directions of Fermi velocity 
vf, V — -vf/vf is the unit vector in the direction of 
vf, — TTT(2n + 1) are Matsubara frequencies, Al^h 
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stands for the order parameter in the left (right) bank, 

-ie + 0) 



Analytic continuation of 501(0) {i.e. io^, 



gives the retarded Green's function; the poles of 5e(0) de- 
termine the energy of the local Andrccv states in the sys- 
tem (i.e. bound states at the interface). The direction- 
dependent Doppler shift pfVs results in the modifica- 
tion of current-phase dependencies and, particularly, in 
the appearance of the countercurrent along the bound- 
ary. The following sections discuss this phenomenon for 
the contacts of two conventional and of two d-wave su- 
perconductors. 



III. WEAK LINK BETWEEN 
CURRENT-CARRYING CONVENTIONAL 
SUPERCONDUCTORS 

In the case of a weak link in the form of a constriction 
between two conventional superconductors with Aj^ = 
Ar = Ao (Ao = Ao (T, V,,)) [7] Eqs. 1-2 can be rewritten 
to study the two coniponc^nts of the current, jx (through 
the contact), and jy (along the contact): 
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where we have split 5^(0) in accordance with the 
two terms in the numerator, which define x- and y- 
components of the current: 



ji(0) = 167reiVot;Frim ^ (^?i5i),,>o,.,>o • (4) 

Wn>0 



Here 12 = 



x,y. We denote jx{0) = jj 



and detach transport supercurrent density jr in jy{0), 
introducing j = jy{0) - jr- 

At T = 0, from Eq.4 after integration we ob- 
tain (see also Fig.l) at ppVs < A^: jj 
jc,osgn{cos |) sin | and j = 0; at < ppVs < Aq: 



jj = jf (^1 - ^ (arccos/3 - , 



(5) 



j=je,osin|(l-/32), /3: 
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where jc,o = f |e| iVoWirAo(r = 0,Vs = 0) is the crit- 
ical Josephson current at T = and Vg = 0; for 
PfVs < Aq the transport supercurrent is linear in v^: 
Jt = ^eNoVFPFVs- The condition ppVs > A^ means 



that (/> G (01, (/)2), where (/>i = 2arccos^2^i '/*2 = 27r— 




FIG. 1. Zero-temperature current-phase dependence for 
Josephson current jj (both in the absence and in the presence 
of the transport supercurrent) and the countercurrent j. The 
value of the transport supercurrent is chosen to be half of the 
critical current (at Vg = 0.5Ao/pf)- 

Thus, at pfVs < A^ the Josephson current is the same 
as in the absence of the transport supercurrent and the 
current density along the contact is equal to the trans- 
port supercurrent density. At pfVs > A^ for a fixed 
value of (j), the Josephson current jj is suppressed by 
the transport supercurrent (compared to the Josephson 
current in the absence of the transport supercurrent): 
jj{vs) < jj{vs = 0), and the countercurrent j appears so 
that the total tangential current density at the contact 
jy(0) consists of the transport supercurrent jr, carried 
by the condensate, and of the countercurrent j, carried 
by nonthermal interface-induced quasiparticlcs. If there 
is a current at the contact which fiows in the direction 
opposite tojhe direction of the current far from the con- 
tact (when j > jx), then the current distribution pattern 
contains vortex- like formations (see also [7]). 

The dependence of the countercurrent j{vs) on Vg 

at T ^ can be deduced from the fact that the current 

increases linearly with Vs for Vs < T/pF and goes to 
zero when Vs tends to the critical value Vg^{T) (^^''(O) ~ 
Pf/Ao). 

Our simple model allows straightforward generaliza- 
tions. For example, consider the situation when the 
boundary transparency is not one, D ^ 1 [11]. Then 

the proximity gap is A^ = Aoy^l — Z)sin^(0/2) and the 

condition of the appearance of the countercurrent (which 
is PfVs > A^ at T = 0) at 7^ 1 is the stronger con- 
dition than the condition at £> = 1. At T = and 
with given Vs and (f> the countercurrent appears for trans- 
parent enough junctions with D > Dc{vs,4>,T = 0) = 



1 



/ sm I . 



Eq.3 determines the energy of Andreev states: £a = 
± A^ — PfVs . Thus there are transport current induced 
zero-energy states which are characterized by the values 
(j) and Vs- At T = the conductance measures the exci- 
tation spectrum of superconductor at the interface [12], 
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[5]; the zero-energy states are known to be responsible 
for the zero-bias conductance peak (ZBCP) [5], [13]. At 
(f) = 0, ZBCP is expected only in the narrow interval of 
the value of the transport supercurrent (when there are 
zero-energy states): Aq < pfVs < 1.03Ao [14]. In the 
system under consideration, where the proximity gap 
is smaller than the bulk value Aq, ZBCP is expected 
even for smaller values of Vg: A^ < ppVs < 1.03Ao. 
The observation of ZBCP can be proposed as a test of 
interface-induced transport-current-dependent quasipar- 
ticles states [5] . Similar conclusions about a conductance 
peak due to the zero-energy bound states were done in 
[15] for the ac Josephson effect in the point contact with 
phase difference close to tt. 

The appearance of the countercurrent can be under- 
stood as the response of the weak link with negative self- 
inductance to the externally injected transport supercur- 
rent. The state of the junction in the absence of the trans- 
port supercurrent at T = is unstable at (/> = tt from the 
point of view that small deviations d(j) = ±0 change jj 
from to Tjc,o- As was shown above, the response of 
the Josephson junction to small transport supercurrent 
at 4> = TT, T = produces the countercurrent j = jc,o- 
It is similar to the equilibrium state with the persistent 
current in ID normal metal ring with strong spin-orbit 
interaction: there is degeneracy at T = and 4> = n, 
and the response of the ring is different at (5(/> / or 
B ^ 0, where B is the effective magnetic field which en- 
ters in the Hamiltonian through the Zeeman term (which 
breaks time-reversal symmetry) [16]. The degeneracy is 
lifted by small effective magnetic field so that the per- 
sistent current rapidly changes from to its maximum 
value. In the case of the weak link between two super- 
conductors in the absence of the transport supercurrent 
there is degeneracy between +py and —py zero-energy 
states; both the time-reversal symmetry breaking by the 
surface (interface) order-parameter and the Dopplcr shift 
(due to the transport supercurrent or magnetic field) lift 
the degeneracy and result in the surface (interface) cur- 
rent [2], [3]. 

IV. WEAK LINK BETWEEN 
CURRENT-CARRYING i?-WAVE 
SUPERCONDUCTORS 

Here we consider the weak link between two d-wave 
superconductors whose axes form an angle 7r/4 [8]. The 
link of this type was proposed as a solid state qubit [10]. 
The equilibrium state of the contact is doubly degenerate 
at 4> = ±7r/2; there is the interface spontaneous current 
js in the equilibrium state; the spontaneous currents flow 
in the opposite directions in the two states: js (— 7r/2) = 
—js One of the ways to control the state of such 

a system is by means of transport supercurrent. When 
the relative angle between two superconductors is equal 



to 7r/4, the current at the contact is described by Eqs.1-2 
with Ai = Aq cos2i? and A/{ = Aq sin2i? (here ^ is the 
azimuthal angle of vp)- The contact can be realized in 
a simply connected geometry, such as a grain boundary, 
or in a doubly-connected geometry such as a ring (when 
the phase difference is controlled by the magnetic flux 
piercing the ring) . Correspondingly we consider here two 
situations: = (no external phase difference) and = 
±7r/2 (in the ground state of the contact). The effect of 
the transport supercurrent on the current parallel to the 
boundary at = is shown at Fig. 2a. 




.0,10 I ■ 1 1 1 ■ 1 1 1 -0 10 I I I I I I I I 

0,0 0,2 0,4 0,6 0,8 q 2 4 6 8 



FIG. 2. (a) tangential current density jy in the banks, 
jT = jyipo), and at the contact plane, jy{Q), at T = O.lTc 
(solid lines) and T = 0.3Tc (dashed lines) versus superfluid ve- 
locity Vs; (b) tangential current density at the interface jy{Q) 
for two values of the phase difference versus Vs (T — O.lTc). 
Here Aoo = Ao(T = 0,v, = 0), jo = 47r |e| N{0)vfTc. 

Making use of Eqs.1-2, we plot the transport supercur- 
rent density jx far from the contact and the tangential 
component of the current density at the contact jy{0). 
The latter flows in the opposite direction to the trans- 
port supercurrent as in the case of contact of conven- 
tional superconductors considered above. For the contact 
of d-wave superconductors, considered in this section, the 
countercurrent appears for close to and tt (in contrast 
to the contact of conventional superconductors, where it 
appears at ^ tt). This countercurrent j = jy{0) — jt, 
which is carried by the quasiparticles, at low tempera- 
ture rapidly grows at small Vg as a function of Vg until 
it becomes equal to the critical value in accordance with 
what was in detail described in [4]. We also note that 
at T = 0.3Tc our results are similar to the results of 
[5], presented there in Fig. 2 (a), where the dependence 
of the tangential current density on the distance from 
the boundary for d-wave superconductor was plotted. At 
this point our conclusions for the existence of the coun- 
tercurrent agree with the previous theoretical results [4] , 
[5] which are consistent with experimental results [1]- [3] . 
However, to the best of our knowledge, there was no sim- 
ilar study for the case of the weak link (particularly, be- 
tween conventional superconductors), controlled by both 
Vs and 0, where the countercurrent along the boundary 
can appear, as described here. 

At = ±7r/2 (see Fig. 2b) the transport supercur- 
rent removes the degeneracy when jy{0) = ijs- The 
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dependence of jy{0) on is nonlinear. The resulting 
tangential current is not always the sum of the spon- 
taneous current js and the transport supercurrent Jt- 
3y{4> = -7^/2) ^ is + jt; iy{(t> = 71-/2) ^ -js+jr for 
Vs < 0.2A/PF, i.e. until \jy \ < j^"'^. 



V. DISTRIBUTION OF THE CURRENT IN THE 
VICINITY OF THE CONTACT 

To illustrate the spatial distribution of the current den- 
sity in the vicinity of the contact we study the case of bal- 
listic point contact between d-wave superconductors (see 
Fig. 3, where thick lines denote the impenetrable parti- 
tion between the superconducting banks). The current 
distribution pattern in the vicinity of the point contact 
between conventional superconductors was considered in 
[7]. 



(a) 

~ ^ ^ J J - * 

^ . 't t .' 


- ! ' r r 


^^^^^ , ^ 


- t \ ^ ^ ^ 
\ V ' - ^ *- 



(hy > * * t ' 

' - ' t 4 li 1 




- - - - t 

^ ^ -* >. ^ 


\ >* ^ -* ^ ^ 
^ ^ " ' ^ 


* " ^ > t * * 


* V ' ' ' ' ' 
\ \ \ i ' ' ' 

1 t » V > ' ' 
t I V * t » ' 



: :| 






' - • I ; / 


\ \ \ v ' - ■ 1 




X W \ > > - T 


' / / ^ / y 


^> \ \ \ \ ^ ' 


' ✓ *^ ^ / ^ 


W ' , 


» w / / / ' 


* f f * I \ >■ , 
It 1/ -i ^ ' 


* >1 ■>! ^ V >. 

^ X w \ \ 


it 1^ ii ^ ' "I 


> >. 1. \ X X 


ii ^ ^ t ' - ■ 1 


- ' li V \ -1, 


i( 1 


- ' ' t V V 




•-■lit 



1 * 1 * ' ' 
; * * i * ' ^ 
+ t i ^ 




% 


\ \ 

k k 
\ \ 


^ 1 1 ( t 
^ * 1 1 1 » t 
» 1 1 1 1 * t 



FIG. 3. Distribution of the current density in the vicinity of 
the contact for (a) g = 0, = 7r/2; (b) g = 0.2 and 4> = 'r/2, 
(c) q = 0.2 and (j) = 3tv/2; (d) q = 0.2 and <f) = tt; here 
q = PFVs /Ao {T = 0,vs=0). 

Although the condition that the contact size is smaller 
than the coherence length a < is hardly realizable for 
high-Tc superconductors, we consider this model as an il- 
lustrative case to show: (a) how the current is distributed 
in the ground state of the contact; (b)-(c) how the trans- 
port supercurrent modifies the current distribution in the 
ground state (qualitatively, the resulting current is a sum 
of what was in the absence of Vs and of the transport 
current); (d) how the appearance of the countercurrent 
results in the vortex-like current distribution. 



VI. CONCLUSION 

The weak link between current-carrying conventional 
and c?-wave superconductors has been studied. In such a 



system it is interesting to study the coexistence of cur- 
rents of different origin, both transport and interference, 
or in other words, to understand how current- carrying 
states of the banks mix non-locally in the vicinity of 
the contact. The system is also interesting because of 
possible applications: in the Josephson transistor with 
controlling parameters (f> and Vg governed by external 
magnetic flux and the transport supercurrent [9], and 
in solid-state qubits, based on a contact of cJ-wave su- 
perconductors [10]. The current at the contact (i.e. its 
components through the contact and along the contact 
plane) is controlled by the values of (p and Vs- For a 
particular relation between (/) and Vg the countercurrent 
along the contact appears. This countercurrent is carried 
by the interface-induced quasiparticles which appear due 
to the coexistence of the transport supercurrent and the 
proximity gap. When the countercurrent density is larger 
than the transport supercurrent density there are vortex- 
like formations in the current distribution in the vicinity 
of the contact. 
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